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The phase structure of Polyakov-loop extended Nambu-Jona-Lasinio 
(PNJL) model is explored at imaginary chemical potential, with partic- 
ular emphasis on the deconfinement transition. We point out that the 
statistical confinement nature of the model naturally leads to characteris- 
tic dependence of the chiral condensate (qq) on 9 — (J,i/T. We introduce 
a dual parameter for the deconfinement transition by making use of this 
dependence. By changing a four-fermion coupling constant, we tune the 
location of the critical endpoint of the deconfinement transition. 

PACS numbers: 11.30.Rd, 12.38.Aw, 12.39.Fe, 25.75.Nq 

1. Introduction 

Phase transitions in QCD have been extensively studied in lattice quan- 
tum chromo dynamics. While recent development enables us to perform 
numerical simulations at physical quark masses, which revealed a crossover 
nature of the QCD phase transition at finite temperature [T], analyses at 
nonzero quark chemical potential fi have been limited to small ^ region 
due to the complex fermion determinant, known as the "sign problem" [2]. 
One of several methods circumventing this problem is to use an imaginary 
chemical potential /j, = i/zj. Indeed, this method has provided transition 
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lines in the T — \x plane via an analytic continuation from those obtained 
at imaginary \x [3113 [5]. Moreover, it has been known that there is a phase 
transition specific to the imaginary chemical potential characterizing the 
deconfinement phase at high temperature [6]. Rich phase structures later 
found in the lattice simulations provide a testing ground for understanding 
the nature of phase transitions in QCD [HJ [TO] . Those properties give 
constraints on model studies which can be extended to real \x. In this work, 
we study the phase structure of the Polyakov-loop extended Nambu-Jona- 
Lasinio (PNJL) model IT5] which satisfies fundamental symmetries of 
QCD relevant for phase transitions at imaginary chemical potential. Focus- 
ing on the deconfinement transition, we show that the "statistical confine- 
ment" feature of the model naturally leads to characteristic behaviors of the 
order parameters while details depend on the choice of the Polyakov loop 
potential. We discuss dual parameters to characterize the phase transitions. 
Finally, we point out the existence of the critical endpoint (CEP) associ- 
ated with the deconfinement transition at imaginary chemical potential and 
clarify the relation between its location and the chiral phase transition. 

In the next section, we will give a brief introduction of the model. We will 
discuss the characteristic behavior of the order parameters as well as the dual 
parameters in Sec. O The critical endpoint of the deconfinement transition 
will be discussed in Sec. H] and Section [5] is devoted to the summary. More 
details can be found in Ref. [13]. 

2. PNJL model at imaginary chemical potential 

The Lagrangian of the two-flavor PNJL model is given by 

£ = q(ij^ - m )q + G s [(qq) 2 + (qi l5 rq) 2 } - U($[A], $*[A}; T). (1) 

The model is an extension of the NJL model, which is an effective model of 
chiral properties of QCD [TOJITO], such that quarks couple with background 
gluonic fields described by a Z(3) symmetric effective potential U which 
takes care of confinement. In the covariant derivative = — iA^, only 
the temporal components of Aq = gA^X 0, /2 is included. The effective poten- 
tial U is expressed in terms of the traced Polyakov loop and its conjugate, 
$ = (Tr c L)/3 and $* = (Tr c L^)/3, respectively. This coupling between 
quarks and gluons leads to an almost simultaneous crossover of the chiral 
and deconfinement transitions at finite temperature, of which order param- 
eters are chiral condensate a = (qq) and the Polyakov loop $ [11] , provided 
the Polyakov loop potential U yields a first order transition at To = 270 
MeV in accordance with pure SU(3) lattice calculations. Two functional 
forms of U, which reproduce the thermodynamic quantities obtained in pure 
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SU (3) lattice gauge theory [12], have been used. One has a polynomial form 

with a set of parameters given in |13j . The other is a logarithmic one [T7] 
= _o^) ^ + ^ _ gM + + _ 3( - M )2j ^ 

The logarithm restricts possible values of <3? and $* to the so-called target 
space, since the argument of the logarithm must be positive. 

At imaginary [/,, the two Polyakov loop variables and <3?* are complex 
conjugate |18j . Moreover, the partition function of the PNJL model at imag- 
inary chemical potential has been shown [18J to have the same periodicity 
in 6 = HijT as that of QCD, Z(6 + 27t/3) = Z(6), which was pointed out by 
Roberge and Weiss [6] as a remnant of Z(3) symmetry. Therefore we may 
express them by using a modulus and a phase = |<I>|e^ and $* = |$|e -1 ^. 

The thermodynamic potential in the mean field approximation reads 

0(T, V, 9) = (G s a 2 + U)V - AV j ^ [3(E p - E° p ) 

+Tln[ 1 +3|$| e ' ; ( e +*)-/ 3E p + 3|$| e i ( 2 e-</>)-2/3£; P + e m-WE P ^ 
+Tln[ 1 +3\$\ e - i ( e+ ^~ l3E p + 3|$| e *M>-20)-2/3£ p + e -3i6-3/3E p ^ / 4 x 



where E p = ^Jp 2 + M 2 , E° = ^p 2 + m%, and M = m - 2G s a. The 
first term in the momentum integral is a divergent vacuum term, which is 
regularized by a three-momentum cutoff A. The cufoff and coupling are fixed 
to G s = 5.498 GeV~ 2 and A = 0.6315 GeV so as to reproduce the vacuum 
pion mass and pion decay constant with mo = 5.5 MeV. In the following 
we mainly focus on the result in the chiral limit too = to preserve the 
chiral symmetry in the Lagrangian. The chiral condensate a serves as an 
order parameter for the chiral phase transition. The order parameters are 
determined by the minimum of the potential which is obtained by solving 
the gap equation dQ / dXi = with Xi = M, \ $\,<p. 



3. Behavior of the order parameters 

3.1. Order parameters at imaginary chemical potential 

First we consider two extreme limits in order to see characteristic 6 de- 
pendences of a which has the same periodicity 2-7r/3 as Expanding Eq. ^ 
for small e~@ Ep , we have a gap equation at small |$| limit in which only a 
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Fig. 1. Left : Chiral condensate for various temperatures as functions of 8. Right : 
Phase of the Polyakov loop 4>. Both results are in the chiral limit tuq = and for 
the logarithmic potential. 
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Fig. 2. Left : Modulus of the Polyakov loop |$|. Right : target space of the 
Polyakov loop on complex $ plane. A region inside the solid lines denote the 
target space of the logarithmic potential in which the argument of the logarithm 
is positive. 

term proportional to cos 39 remains with a small magnitude ~ e ~ 3 l 3E p indi- 
cating the statistical confinement. This dependence naturally leads to the 
periodicity 2tt/3 when |<I>| is negligible and chiral symmetry is broken. On 
the other hand, when |<3?| ~ 1, the model reduces to the NJL model except 
for the coupling of </> with 9, as seen in Eq. In this case, the apparent 
9 dependence is governed by cos 9 as a consequence of deconfinement. Al- 
though this factor does not match with the required periodicity 2tt/3, it is 
preserved by a change of 0, namely, the Roberge- Weiss transition. 

We show the chiral condensate in the left panel of Fig. [1] obtained by 
numerically solving the gap equations. One sees that a at low temperature 
(T = 220 MeV) exhibits small and smooth variation as a function of 9, as 
discussed above. On the other hand, one sees a cusp at 9 = ir/3 and T = 280 
MeV. This is a consequence of a Roberge- Weiss transition depicted in the 
right of Fig. [TJ in which the phase eft changes from to — 2n/3, smoothly at 



morita printed on November 21, 2011 



■5 



Log, m =0 




Chiral 




■ CT=0 






RW 


<j - cose 


- cos(2it/3-6) 


Deconfinement 




o 


- cos39 



0.34 

5" 0.32 
a) 

£ 0.3 
I- 

2. 0.28 

I 0.26 
a. 

§ 0.24 
I- 

0.22 
0.2 



0.5 



1 

9/(it/3) 



1.5 



■ Polynomial, m =0 




Chiral-' 

■ a=0 


a=0 


"" a~cos8-"' 


RW *•-. 
cT*-co5(2jt/3-e")" "■ 


Deconfinement 






- cos39 



0.5 



1 

e/(it/3) 



1.5 



Fig. 3. Phase diagram on T — 9 plane. Solid lines, dotted lines, and dashed lines 
stand for first, second, and crossover transitions, respectively. 



low T but discontinuously at high T. As a result, the required periodicity 
of a is preserved. One also sees a second order chiral phase transition for 
T = 280 MeV in which the chiral symmetry is broken around 9 = ir/3. This 
implies the chiral critical temperature at imaginary fi is higher than that of 
zero and real [i. This can be also understood from the gap equation for a, 
since cos n9 is replaced by cosh nf3/j, for real ii. 

While the above properties are independnent of the choice of U, there 
are some potential dependent features as follows. In Fig. Q] and the left 
of Fig. [21 one sees a discontinuity in the order parameters at the same 9. 
This shows a first order deconfinement transition which exists only in the 
case of the logarithmic potential ([3]). The polynomial potential exhibits 
smoother change near phase transition. In the right of Fig. [21 the target 
space of the Polyakov loop is displayed. Owing to the Z{2>) symmetry, U has 
three degenerate minima at T > T$. Putting quarks into the system makes 
one of those minima favored. While Im<I> = is always chosen at 9 = 0, 
Im<I> is favored at imaginary chemical potential due to the coupling 
of 9 and <j) seen in Eq. At low temperature where minimum of U is 
close to the origin, the minimum of the effective potential smoothly moves 
from (f) = to 4> = — 27r/3 across 9 = tt/3. At high temperature, however, 
there is a potential wall which makes the transition from point A to point B 
discontinuous. Since the polynomial potential does not have any restriction 
of the target space, the minimum passes outside (C) the target space near 
the RW transition. 

The phase diagrams shown in Fig. [3] summarize the behavior of the 
order parameters. One sees a first order deconfinement transition and an 
associated critical endpoint (CEP) only for the logarithmic potential. This 
also implies that the RW endpoint, where the first order RW transition 
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Fig. 4. Dual parameters compared with Polyakov loop. Top panels show <I>, E^ 1 - 1 
and Eg for U\ og (left) and U po \ y (right) while bottom ones displays their derivatives 
with respect to temperature. 



terminates, is a triple point. On the other hand, one sees a second order RW 
endpoint for the polynomial potential. The properties of the RW endpoint 
in QCD might reflect the nature of the QCD phase transition at real [/, We 
refer to Refs. [8] and [7] for recent calculations of Nf = 2 and Nf = 3 latttice 
QCD, respectively. Especially it should be noted that the order of the RW 
endpoint has a non-trivial bare quark mass dependence which cannot be 
reproduced by chiral effective models. (See Sec. S|) An improved model was 
proposed in Ref. [19J to reproduce this property. 



3.2. Dual parameters for deconfinement 



It has been shown that information on the deconfinement is encoded in 
9 dependence of a. We can consider dual parameters which characterize 
the deconfinement transition. A dual parameter was introduced in [20] . By 
considering a twisted boundary condition for quarks q(x, f3) = e lv q(x,0), 
one may define the corresponding chiral condensate cr(<p). Then the dual 
chiral condensate T,( n ' reads 



2~ 



dip 
— < 
2vr 



-imp 



i (Tr[(m + D v y 



(5) 
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While the twisted boundary condition is similar to introducing imaginary 
chemical potential [21], it does not apply to the background gauge field. 
Therefore, cr((p) has a periodicity 2n and was calculated in a PNJL model 
by fixing the Polyakov loop at its 9 = value [22J. Particularly E^ 1 ) is 
called dressed Polyakov loop, since it has the same transformation prop- 
erties under Z{3) and thus is expected to serve as an order parameter of 
the deconfinement transition. Analogously, we consider a modified dual 
parameter which utilizes the characteristic property of cr(0), 

4 n) ( T ) = T- d6e- md a(T,e). (6) 

2VT J-n/3 

where we take the integration range [— 7r/3, vr/3], owing to the periodicity 
of <x(0). 

We compare those dual parameters for n = 1 with the Polyakov loop 
in Fig. HI as well as their derivatives with respect to temperature, of which 
peaks can be regarded as (pseudo) critical temperatures. One sees that while 
dual parameters show a rapid increase as seen in the Polyakov loop (top) J]] 
their derivatives exhibit different peak structures. The derivatives of the 
dual parameters have a peak at the chiral transition temperature, inde- 
pendent of U. As for the deconfinement, however, existence of the peak 
depends on U. The dressed Polyakov loop exhibits a peak for the U\ og for 
which |<3?| shows stronger crossover than lA vo \ y . Moreover, the modified dual 
parameter exhibits only a shoulder even for U\ og . This result indicates dif- 
ferent sensitivity of the dual parameters to the chiral and deconfinement 
transition. 



4. Critical endpoint of deconfinement 

Now let us turn to the deconfinement CEP found in the case of U\ og . 
Here we vary the four-fermion coupling constant G s to preserve the chiral 
symmetry in the Lagrangian. Locations of the CEP are shown in the left 
of Fig. [5] for various values of G s . One sees the squared critical chemical 
potential ^cep increases with G s to reach ^ EP = around G s ~ 6.3 
GeV -2 . In the right of Fig. [5j we also depict a phase diagram for G s = 6.5 
GeV -2 in which the CEP exists at real chemical potential. One sees that 
the first order deconfinement transition starting from the RW endpoint (see 
Fig. [3]) is prolonged, while the chiral critical line moves upward. The relation 
between these two changes can be understood as follows. Since the Polyakov 
loop potential U\ og has a first order phase transition at T = Tq, the model 
results in the same transition when the effects of quarks are negligible in 



1 Dual parameters are normalized to as T — > and 1 as T — > oo |14] , 
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Fig. 5. Left : location of the deconfinement CEP (Tcep, A*cep) as a function of G s 
Right : Phase diagram for G s = 6.5 GeV~ 2 which gives ^cep > 0. 
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Fig. 6. Behavior of |$| for large G s . Left : U\ og . Right : U po \ y 



thermodynamics. The contribution of quarks to thermodynamic potential is 
essentially determined by the dynamical quark mass M = m,Q — 2G s a, not by 
the current quark mass mo, as seen in Eq. When dynamical quark mass 
becomes lighter around T = Tq, the deconfinement transition is modified 
to a crossover one. As G s increases, the stronger coupling leads to a larger 
condensate \cr(T = 0)| thus the dynamical quark mass becomes heavier. 
This appears as the modified chiral critical line in the phase diagram at 
G s = 6.5 GeV -2 and the resultant dynamical quark mass is heavy enough to 
recover the first order the deconfinement transition. At the reference value 
of G s , the imaginary chemical potential weakens the thermal terms by cos n0 
in the thermodynamic potential thus resembling a heavier quark mass which 
yields the CEP and a first order transition. While the above consideration 
is completely independent of the form of U, quantitative features such as 
the value of dynamical quark mass which makes the transition first order 
depend on the choice of hi. 

Figure [6] shows the behavior of |4>| for various G s at vanishing chemical 
potential. One sees that |$| becomes steeper for larger G s in both of U. The 
case of U\ g has a discontinuity already at G s = 6.5 GeV -2 as mentioned 
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above. U po \, which has a smoother variation of |$| against T, eventually 
approaches the pure gauge case for much larger G s . At G s = 25 GeV -2 , 
where the dynamical quark mass at T = is around 2.5 GeV, the first 
order deconfinement transition is recovered. The origin of this difference 
is the much weaker first order transition in U po \ y , which easily turns into 
crossover when quarks heavier than 2.5 GeV are put into the system. If 
one characterizes a strength of the deconfinement transition by a gap of 
the Polyakov loop A<£> at T = T , one finds A$ = 0.47 for U iog and 0.072 
for U po \j. Since G s determines the scale of the dynamical chiral symmetry 
breaking, a(T = 0), this result indicates an interplay of the two transitions 
which have a unique scale, Aqcd, in the case of QCD. 

It has been shown that a first order deconfinement phase transition also 
emerges in the large N c limit of the PNJL model [23]. This result has a 
common origin with the present study in a sense that taking large N c limit 
makes the system gluon dominated due to 1/N C suppression of the quark 
contribution, while large G s thermally suppresses quarks in the chirally 
broken phase. In our case, however, chiral transition temperature moves 
upward thus there is a discrepancy between the deconfinement and chiral 
transition temperatures, in contrast to the large N c limit with a fixed G S N C 
in [23] . 

5. Summary 

We have explored the deconfinement transition in the PNJL model at 
imaginary chemical potential. We point out that the chiral condensate at 
imaginary chemical potential, c(9), has a characteristic 9 dependence due to 
the deconfinement property which naturally arises from the statistical con- 
finement feature of the model. While the confined phase is characterized 
by a smooth cos 39 dependence, the deconfined phase exhibits cos 9 depen- 
dence together with cusps at 9 = ir/3 (mod 27r/3) induced by the abrupt 
change of the phase of the Polyakov loop (Roberge- Weiss transition). We 
introduce a new dual parameter utilizing this 9 dependence and compare it 
with the Polyakov loop and the dressed Polyakov loop. Different sensitivi- 
ties of these parameters to chiral and deconfinement transitions are found. 
Changing the four fermion coupling constant, we found that an interplay 
between the thermal quark contribution through the dynamical chiral sym- 
metry breaking and the Polyakov loop potential determines the location 
of the deconfinement CEP at imaginary chemical potential. In particular, 
we found that the deconfinement CEP can be located in the real chemical 
potential regime for a Polyakov loop potential with a strong first order tran- 
sition and large dynamical chiral symmetry breaking. We expect that these 
results are useful for understanding of the QCD phase transition. 
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